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Abstract 

O ; 

^ A short review of the optical and far-infrared measurements on MgB2 is given. Multiband and multigap effects are analyzed 
Q by comparing optical properties with other experiments and ab initio calculations. The covered topics are: the plasma frequency, 
^ electron-phonon interaction, impurity scattering, the effects of C and Al substitution, interband transitions and the far-infrared 
C/3 signatures of the superconducting gaps. 
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1. Introduction 

The existence of two distinct superconducting gaps in a 
combination with a comfortably high T c of 40 K in MgB2 
offers a unique opportunity to study the electrodynamics 
of a multigap superconductor. The effects of interband cou- 
pling and scattering on the superconducting order param- 
eter and T c were envisaged long time ago [TJ , but only after 
the discovery of superconductivity in magnesium diboride 
[5] the community was given a real chance to study the phe- 
nomenon experimentally. MgB2 seems to be a very unusual 
case of 'conventional' superconductivity, where the novel 
multiband and multigap physics can be well captured by ab 
initio calculations of the electron and phonon dispersions, 
electron-phonon interaction and properly modified strong 
coupling Eliashberg formalism [3]. As a result, a rather di- 
rect quantitative comparison between theory and various 
experimental probes is plausible, in contrast to the case of 
high-T c cuprates with an important role of electron corre- 
lations and still debated mechanism of superconductivity. 
' Electronic band structure is the starting point to 
understand the unusual superconducting properties of 
MgB2 . The first- principles calculations that exist since the 
late 1970's [4] were refined after 2001 by several groups 
|5|6|7|8|9j . Angle-resolved photoemission (ARPES) [TO] 
and de Haas- van Alphen (dHvA) measurement [11] showed 
a good agreement with the calculated band dispersion. 
The two-gap behavior, first predicted theoretically [T2"] . 
was confirmed by tunneling spectroscopy [113 14 , specific 
heat [TS], ARPES [IB], Raman [T7] and other techniques. 

Until recently, the agreement of optical data with the- 



oretical predictions has been far less successful. Perhaps, 
the most notable and serious issue was an extremely small 
plasma frequency reported from several early optical stud- 
ies |18I19I20I21I2"2] . More recent studies [23124125] showed 
that the plasma frequency is actually much closer to the ex- 
pected value, although optics results call for certain amend- 
ments to the theory. The main problem in obtaining rea- 
sonably reproducible optical data is the very small size of 
the existing high-quality single crystals [26127] . 

The electron-phonon interaction, which is responsible for 
the superconductivity in MgB 2 as suggested by the isotope 
effect |28I29] , inevitably manifests itself via the energy de- 
pendent renormalization of the scattering rate and the ef- 
fective mass. The latter parameters can be derived from 
the optical spectra already in the normal state, provided 
that the correct plasma frequency is available. Single crys- 
tal measurements [23 have shown a reasonably good agree- 
ment with the ab initio calculations |6I12I9| and confirmed 
a stronger electron-phonon coupling in the a bands. 

While the biggest controversies regarding the general 
normal-state optical properties of MgB 2 seem to be settled, 
the main issue related to far-infrared spectra in the super- 
conducting state, namely the absence of a clear evidence 
of the two-gap structure, still remains. Unlike the situation 
with early optical measurements on magnesium diboride, 
a number of different far-infrared experiments, especially 
those on thin films, are quite consistent with each other. 
Therefore we have to ask ourselves whether the existing 
theory describes adequately the infrared optical response 
of a two-gap superconductor. 

The organization of this paper is as follows: In Section 
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[5] some basic considerations regarding the multiband con- 
ductivity are given. Section [3] deals with the experimental 
determination of the plasma frequency. In Section [4] the 
electron- phonon and impurity scattering as observed by op- 
tics are discussed. Interband transitions are considered in 
Section [5l Finally, the far-infrared signatures of the super- 
conducting gap(s) are reviewed in Section [SI 

2. Conductivity of a multiband system 

The two main contributions to the complex optical con- 
ductivity cr(u>) = a\{uj) + ioiiui) of a metal are the intra- 
band (Drude) component centered at zero frequency and 
the one due to momentum conserving interband transitions 

a(u>) = ctd(cj) + cr/ B (o;). (1) 

The detailed theory of multiband transport based on the 
Boltzmann equation was developed in Ref.[30]. A widely 
used approximation for the Drude conductivity <td(u;) of 
a system with several conduction bands is the parallel- 
resistors formula, where each resistor (or better say conduc- 
tor) corresponds to a separate band. Although this seems 
like assuming that quasiparticles in different bands do not 
'see' each other, the impurity or phonon induced interband 
scattering probabilities to a certain approximation can be 
just added to the intraband ones, resulting in a parallel 
band conduction with effectively renormalized scattering 
rates. 

The metallic properties of MgB2 are determined by two 
distinct types of electronic bands coming almost entirely 
from the boron states: the strongly covalent almost two- 
dimensional a bands formed by hybridized sp x p y orbitals 
and the more isotropic it bands made of p z orbitals. The 
holes in the a bands are strongly coupled to the in-plane i?2 ff 
phonon mode at ~ 75 meV, giving rise to a high clcctron- 
phonon coupling constant |6|12|9j . In this paper we shall 
make no distinction between bands of the same type, rep- 
resenting the conductivity by the sum of the two terms 

a D (uj) = ^ ^A/st^O- ( 2 ) 
The optical sum rule reads as follows 

oo 

a 1D {uj)duj = - ^,(3 = g w p: ( 3 ) 

0=<r,n 

where uj p ^ and u> Pj7r are the unscreened (bare) plasma fre- 
quencies of the a and tt bands and uj p is the total plasma 
frequency Of course, the conductivity and the plasma fre- 
quency are both tensors with two different components: 
one parallel to the ab plane and another along the c axis. 
The first-principle calculations predict the a band plasma 
frequency to be very small along the c axis. 

It follows from Eq.([2]) that there is no rigorous exper- 
imental way to separate aD,a{u) and crc i7r (cj). However, 



the multiband structure may become apparent in optical 
spectra due to the strong disparity between the two bands. 
Such a contrast can stem from (i) different anisotropies of 
the plasma frequencies u> p , a and to p ^, (ii) manifestly dif- 
ferent electron scattering on phonons and impurities, and 
(iii) a large difference of superconducting gap values. In the 
following sections these possibilities will be analyzed in de- 
tails. 

3. Plasma frequency 

The plasma frequency that can be measured using the 
optical sum rule ([3]) is determined by the electronic band 
dispersion Ckp [12] : 




(a = a,c) (4) 



V is the sample volume, / is the Fermi-Dirac distribution 
function and f3 is the band index. At not too high temper- 
atures, only states close to the Fermi energy Ep contribute 
to Eq.®, and iv p ~ N(E F )(v F ), where N(e) is the density 
of states and vf is the Fermi velocity. 

The plasma frequency is an important parameter for test- 
ing the consistency of band structure calculations with op- 
tics. Therefore, the reported value of only 1.5 - 3.0 eV for 
ui p by several groups [18 19 20 21 22 soon after the discov- 
ery of superconductivity in MgB2, as compared to the the- 
oretical prediction of 7 cV in both directions, was rather 
puzzling. The situation changed recently, when three in- 
dependent measurements on high-quality single crystals 
[23 24 25 showed the values of cu p much closer to the cal- 
culated ones. 

Given the large spread of optical data, one should crit- 
ically consider the experimental issues. Possible complica- 
tions are (i) sample purity, (ii) surface contamination, (iii) 
analysis ambiguity due to a limited spectral range and (iv) 
an overlap between the Drude peak and interband peaks. 

Impurities in small amounts are not expected to affect 
significantly the plasma frequency. However, they increase 
the elastic scattering and shift the Drude spectral weight 
to higher frequencies. Interestingly, in MgB2 impurities can 
affect differently the charge scattering in a and n bands 

As observed by several groups, the optical quality of the 
surface of MgB 2 poly- and single crystals degrades in air 
within minutes after polishing |19|21|23|25j . To the author's 
knowledge, there is no information about the surface degra- 
dation of MgB2 films. The surface contamination strongly 
reduces the absolute reflectivity, pushing down the appar- 
ent value of the plasma frequency. 

A source of inaccuracy not to be neglected is a set of 
assumptions taken to experimentally derive and integrate 
tTi(w). In particular, the error bars due to high- and low- 
frequency extrapolations in the Kramers-Kronig (KK) 
analysis of reflectivity that are hard to estimate may af- 
fect significantly the value of u> p , especially if spectra are 
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available only up to few electronvolts. This uncertainty is 
reduced if the ellipsometric technique is used allowing the 
direct measurement of <j\ (lo) and oi (lo) . 

Tu et al. [18 measured the reflectivity of high-quality 
oriented films (T c =39.6 K) with the c axis perpendicular 
to the surface in the range from 3.5 meV to 2.7 eV and 
extracted <j(lo) by the KK transform. A value of ui p ~ 1.8 eV 
was deduced from the sum-rule analysis. Notably, the film 
reflectivity is very high at low frequencies (compatible with 
the predicted high metallicity of MgB2) but becomes quite 
low above 0.5-0.7 eV as compared to the recent data on 
single crystals 23 24 25] . The same trend of high-frequency 
reflectivity suppression is present in the film study by Mun 
et al. [20] and, to a larger extent, by Chvostova et al. [22 . 
It is not clear at the moment whether this is due to surface 
contamination, film strains or other reasons. 
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Fig. 1. Optical anisotropic spectra of MgB2 at 300 K derived from 
ellipsometry and reflectivity measurements on an ac-oriented single 
crystal: the normal-incidence reflectivity R(u)(&), optical conductiv- 
ity (Ti (oj), the dielectric function ei(ui) (b) the effective number of 
carriers N e g (uj) and the effective plasma frequency f2 p ff((^) as a 
function of the integration cutoff energy (c). Adapted from Guritanu 
et al. [23]. 

In another early study, Kuzmenko et al. 19] obtained a 
direction-averaged conductivity of a dense polycrystalline 
sample (T c =39 K) from 6 meV to 4.6 eV by a combination 



of reflectivity and ellipsometry measurements. The conduc- 
tivity shows a narrow Drude peak with a plasma frequency 
of only 1.4 eV and additionally a broad infrared contin- 
uum with a spectral weight corresponding to the plasma 
frequency of about 5 eV. Since this sample has a sizeable 
contamination of MgO, a possible explanation of the ap- 
parently small Drude plasma frequency is that the narrow 
peak is due to the a bands, which are not strongly affected 
by impurities, but renormalized by the electron-phonon in- 
teraction, while the broad continuum is largely formed by 
the 7r bands [32] • This is probably the main reason for the 
low reflectivity of that sample, although some influence of 
the surface contamination cannot be excluded. 

The first attempt to distinguish oo Pta and lo p _ c was un- 
dertaken by Fudamoto and Lee 21j. From a comparison of 
the reflectivity spectra measured on a mosaic of ab-plane 
oriented crystals and on a polycrystalline sample they cor- 
rectly deduced that the plasma edge parallel to the ab plane 
is at about 2 eV while the one along the c axis is at 2.75 
eV. However, the conclusion of Ref. [21] about a significant 
difference between Lo p ^ a and u> PjC was not confirmed by the 
latest single-crystal studies [23125] . The reason is that the 
position of the plasma edge corresponds to the screened 
plasma frequency lo* ~ lo p / ^/£oo which depends not only 
on Lo p but also on the dielectric constant eoo from inter- 
band transitions. The latter quantity appears to be strongly 
anisotropic in MgB2 (Section [5]). 

A reliable determination of the optical anisotropy re- 
quires single crystals with the c axis parallel to the sur- 
face. The major difficulties are the small dimensions of the 
existing single crystals, especially along the c axis (< 200 
/jm) and the mentioned surface degradation. Guritanu et 
al. [23] succeeded in doing spectroscopic ellipsometry in the 
range 0.75 - 3.7 eV and reflectivity from 0.1 to 0.85 eV on 
the ac- and ab-surfaces of small single crystals (T c = 38 
K) at room temperature (Fig[T]), by adding focusing lenses 
to a standard ellipsometry setup and keeping samples in 
a nitrogen flow. Both in-plane and c-axis spectra exhibit 
a metallic behavior, characterized by a reflectivity plasma 
edge, a Drude peak in u\ (lo) and a negative e\ (lo) . However, 
one can see a strong anisotropy of the plasma edge, con- 
sistent with the assignment made in Ref .[21]. The plasma 
frequency derived by the sum rule is about 6.3 eV for both 
directions, if the integration is continued up to 1.7 eV. This 
cutoff value is high enough to comprise most of the intra- 
band spectral weight, but it is below the interband peak at 
2.6 eV along the a axis (discussed in Section[5|). Since ellip- 
sometry provides directly both real and imaginary parts of 
the dielectric function, the value of the plasma frequency 
is practically independent on extrapolations. The value of 
6.3 eV is much closer to 7 eV predicted by the LDA calcu- 
lations. A very small anisotropy of the plasma frequency is 
also in a good agreement with the theoretical predictions. 

Independently, Kakeshita et al. [25 measured a- and c- 
axis reflectivities of MgB2 from 75 meV to 3 eV on ac- and 
ab-oriented single crystals at room temperature using a mi- 
croscope spectroscopy technique (Fig[2]). The sample was 
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Fig. 2. The reflectivity spectra of Mg^i—zd^ for E \\ ab plane (a) 
and E \\ c axis (b). Inset: the x dependence of the plasma frequencies 
LU Pi a and L0 Pt c, and the scattering rates -y a and 7 C obtained by a 
Drude-Lorentz fit. Adapted from Kakeshita et al. [25) . 
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Fig. 3. The optical conductivity spectra [(a), (c)] and the loss function 
[(b), (d)] of Mg(Bi_ 2 ,C a; )2 for E ||ab plane and c axis obtained by 
the KK analysis of reflectivity curves from Fig|2] Adapted from 
Kakeshita et al. I25I . 

polished prior to each measurement in order to avoid the 
effect of the surface degradation. The data from Refs. 25J 
and 23J match rather well with each other, although the 
plasma edges are sharper in Ref . |25j , indicating a narrower 
Drude peak (Fig(3J). The plasma frequencies estimated in 
Ref. [25] from the optical sum rule with a cutoff at 2.2 eV are 
u p ,a ~ 8.4 eV and w PjC ~ 7.0 eV. However, for such a high 



cutoff the interband transition at 2.6 eV contributes signif- 
icantly to the sum rule along the a axis. A Drude-Lorentz 
fit that treats the interband peak as a separate contribu- 
tion gives a smaller value ui p , a ~ 7.4 eV. The values of the 
plasma frequencies obtained in Refs. [53] and [5S] seem to 
agree each other within the error bars of the KK analysis. In 
both cases, oj p ^ a might be slightly overestimated since the 
theory predicts the existence of a weak a — > a interband 
transition at 0.3-0.5 eV. Within the error margins, there is 
no prohibitive disagreement between optics and LDA cal- 
culations, as has been shown by independent measurements 
on single crystals grown in different groups. 
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Fig. 4. Anisotropic plasma frequencies of a and tt bands as a function 
of the Fermi level in a rigid-band approximation. Also shown are 



the total plasma frequencies ui Pt , 
Ref. [3]. 
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The influence of Al, C and other substitutions on the 
properties of MgB2 has been recently the focus of intense 
research. All substitutions are found to suppress T c (|33J 
and references therein). It was argued [33] that this sup- 
pression can be well described by a combination of the fill- 
ing of a bands and interband scattering. Optics contributes 
to a better understanding of the substitution effects. While 
Al replaces Mg and carbon substitutes boron atoms, both 
of them dope electrons to the system, which is expected 
to affect the plasma frequency. In the rigid band picture, 
where the doping merely shifts the Fermi level, the change 
of the band plasma frequencies can be derived from the 
band structure calculations for the pure compound. Figf4] 
shows the plasma frequencies of a and 7r bands as well as 
the total ones for the two polarizations 3J. With the in- 
crease of the Fermi level, both uj p ^ a and o-> PiC are expected 
to diminish, but the decrease of Lo Pya is faster due to a rapid 
decrease of uj p ^ aa - The latter effect is caused by a fast de- 
pletion of the density of states (DOS) of a bands. 

Kakeshita et al. [5S] measured reflectivity on a series of 
carbon doped samples Mg(Bi„ x C x )2, < x < 0.1 (Figf2]). 
Remarkably, the c-axis plasma edge as well as the corre- 
sponding peak of the loss function (Fig[3j) shift to low fre- 
quencies with doping while the one for the a axis is al- 
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most doping independent. Correspondingly, u> PtC notably 
goes down, while w p , a shows only a minor decrease with 
C substitution (shown as inset of Fig[2]). The decrease of 
both ujp^ a and w PiC with doping qualitatively agrees with 
their calculated Fermi-level dependence (Figj4j) . However, 
the rigid band approximation does not explain that u PiC 
goes down much faster than u> Pja . It is necessary therefore 
to take into account the doping induced band modification. 
The calculations of Refs. [33] and [34] indeed show a much 
slower decrease of the tr-band DOS compared to the rigid 
band model. This results in a slower rate of the suppression 
of T c with doping, in agreement with observations. 
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Fig. 5. In-plane reflectivity spectra of C- and Al-doped single crystals 
of MgB2. Adapted from Di Castro et al. [24] . 

The ab-plane optical reflectivity of a series of C- and Al- 
doped single crystals of MgB2 was measured by Di Castro 
et al. 24J (Fig[5]) . The data indicate that each of these sub- 
stitutions causes only a minor shift of the a-axis plasma 
edge. Correspondingly, a quantitative analysis showed a 
slight decrease of the total in-plane plasma frequency. As 
far as carbon doping is concerned, this agrees well with the 
results of Ref . [55] . In Ref. [23] a value of Lo p ^ a of about 5 eV 
was reported, which is somewhat smaller than the values 
obtained in Refs. [53] and [55J. The low- frequency reflectiv- 
ity is close to the one from [53] , but drops faster above 2 eV, 
which partially explains, together with uncertainties due to 
a limited spectral range, the difference of the plasma fre- 
quencies. Another reason is that the model dielectric func- 
tion used by the authors [24] to fit reflectivity contained 
an extra mid-infrared (MIR) peak ascribed to the a — ► a 
interband transition. 

4. Electron-phonon and impurity scattering 

The information about electron-phonon interaction and 
impurity scattering can be derived from the shape of the 
Drude peak. In particular, the strength of the electron- 
phonon coupling, which is widely believed to be responsi- 
ble for the superconducting pairing, can be deduced from 



the infrared spectra. The knowledge of impurity scattering 
rates can be used to test the hypothesis about a strong dis- 
parity between a and it bands and a very small interband 
scattering used to explain the survival of two distinct su- 
perconducting gaps [32] . 

The broadening is conveniently described via the 
memory function M(oj,T) = l/r(w,T) — iu>\(u>,T), 
where 1/t(lu,T) is the scattering rate and X(uj,T) = 
m*(uj, T)/rrib — 1 is the electron mass renormalization (nib 
is the bare band mass). Note that the memory function 
definitions present in the literature differ by a factor of i 
[35 3(3]; we follow the notation, where l/r(w,T) is given 
by the real part of M(ui,T). In the case of a single band, 
M{lo,T) can be directly extracted from the intraband 
optical conductivity via the extended Drude formalism 



M(lo,T) = 



1 



47T a D (u>,T) 



(5) 



For two bands, the situation is more complicated. The 
band-dependent scattering has to be accounted for by tak- 
ing separate memory functions AI a (uj, T) and M^w, T) 



o\d. 



■<»> = h T, 



An ' Mr(lo. T) — iuj ' 

0= CT,7T ^ 



(a = a, c). 



(6) 



In fact, this complication is beneficial since the scattering 
dissimilarity distinguishes optically the contributions from 
different bands. 

We assume that impurities/vacancies and phonons con- 
tribute additively to the memory function (Matthiessen's 
rule): 



Mp(oj,T) = jp timp + Mf3. ph (uj,T), (P = cr,7r). 



(7) 



The impurity scattering is strongly sample dependent and 
frequency and temperature independent, in contrast to the 
one due to lattice vibrations. The expression for the phonon 
memory function at finite temperatures was obtained by 
Shulga et al. [37] 



Mp tPb (u;,T) 



-2i J a 2 Mn)K 
'o 



2ttT' 2ttT 



n \ 



dQ, 



(8) 



which is a generalization of the famous T = formula due 
to Allen [58] . Here af r Ff3(U) is the transport Eliashberg 
function, and 



K(x, y) = — h [Wl — ix + iy) — -0(1 + iy)] + 

y x 



y + x 



[0(1 - ix - iy) - ^(1 - iy)] 



where ip(x) is the digamma function. The total electron- 
phonon coupling strengths are given by the constants 



Kt,b — 2 



(9) 
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Eq. ((SJ) follows from the Kubo formula in the weak-coupling 
limit, although it remains accurate even at A tr ~ 1 [39 . 
There are several ab initio calculations of the Eliashberg 
functions |6|9|12j which all show a strong coupling in the a 
bands and a moderate coupling in the 7r bands. 

The memory function as determined from the experi- 
mental data by Eq.([5]) is proportional to the plasma fre- 
quency squared. Therefore the strong underestimation of 
uip in the first optical experiments resulted in anomalously 
small electron-phonon coupling constants [40] . For exam- 
ple, in Ref. [18] a value A ~ 0.13 was found, which is in- 
compatible with the T c of 40 K without invoking exotic 
coupling mechanisms |41j . 
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Fig. 6. Extended Drude analysis of the optical conductivity of MgB2 
at 300 K along the in-plane and the c-axis directions. The symbols 
are the data, the solid curves show a two-component calculation. 
Adapted from Guritanu et al. |23| . 

FigJSJshows 1/t(w) and m*(w)/mj at 300 K obtained by 
Guritanu et al. [23] by applying Eq. ([5]) to the single crys- 
tal data using cj p a w w P)C ss 6.3 eV. Since the plasma fre- 
quency of cylindrical a bands is predicted by the LDA cal- 
culations [5142] to be very small along the c axis (FigH|), 
the extended Drude parameters along the c axis is almost 
entirely determined by the ir bands. In contrast, the contri- 
butions from a and n bands to the in-plane properties are 
comparable, therefore both scattering rate and mass renor- 
malization for E || ab are effectively averaged over the two 
bands. This is why a larger m* (to) / mb and a faster growth 
of l/r(w) with photon energy for the in-plane polarization 
indicate, without any modeling, that the electron-phonon 
coupling in the a bands is stronger than the one in the tt 
bands. One has to keep in mind though that the a — > a 
interband peak is predicted to be at ~ 0.3-0.5 eV for the in- 
plane polarization (see Section [5]) that may interfere with 
the extended Drude analysis. 

Kakeshita et al. [35] found that the c-axis optical con- 
ductivity is much better described by a single Drude peak 
with a frequency independent scattering rate than the con- 



ductivity parallel to the ab plane. This agrees well with 
the much weaker frequency dependence of l/r(w) for E || c 
than for E \\ ab [23] (FigE|). 

A more quantitative comparison between experiment 
and theory can be made using Eqs.©, (0 and ([8|). Al- 
though not all parameters that enter these formulas (four 
plasma frequencies, two impurity scattering rates, and two 
Eliashberg functions) can be independently derived from 
the data fit, the measured anisotropy of optical spectra 
imposes tight constraints on them. At frequencies above 
the phonon range the data are sensitive to the total cou- 
pling constants but not to the spectral structure of the 
Eliashberg functions. It was found 23J that the ab initio 
calculated coupling constants (~ 1.1 for the u-bands and 
~ 0.55 for the 7r-bands [5|6] ) agree well with the spectral 
data, if the quasi-2D nature of the a bands is assumed. 

It appears that both the in-plane and c-axis infrared re- 
flectivity spectra are quite sensitive to the 7r-band impurity 
scattering , allowing an accurate determination of r y n< imp 
(7^,im P ~ 85 meV for the single crystal from Ref. [23]). In 
contrast, 7 CTi i mp cannot be reliably derived from infrared 
spectra alone, because the scattering in the a band is dom- 
inated by the electron-phonon interaction. The ambigu- 
ity can be fixed by considering the temperature depen- 
dence of the in-plane resistivity p a (T), for which the ra- 
tio p a (300K) I 'p a (40K) (RRR) is rather sensitive to -f a ,im P - 
In order to match the experimental curve measured on a 
sample grown in the same conditions that showes RRR « 
5, 7<T,im P is about 7 times smaller than 7 T .im P had to be 
assumed. Such a larger difference of the band scattering 
rates agrees with Raman data [31] . One should note that 
this supports the hypothesis of Mazin et al. [32] that the 
a holes with wavefunctions strongly confined within the 
boron planes should be scattered by typical lattice defects 
(except the carbon atoms) much less than the 7r-band carri- 
ers. This is due a larger probability of the defect formation 
in the interstitials than inside the covalently bonded boron 
sheets. 

Kakeshita et al. [35] observed a strong increase of the 
scattering rate in the a bands and a significant but lesser 
increase in the n bands with carbon doping. This can be 
explained by the fact that C substitutes boron atoms, af- 
fecting strongly the a band scattering. A similar result was 
obtained by Di Castro et al. [53], who found that carbon 
doping has a more pronounced effect on the a band scat- 
tering than the Al substitution. 

5. Interband transitions 

The interband optical conductivity is determined by the 
band dispersion and transition matrix elements [42)43] 

ffuB ' o(w) = ^^ £ \(kf\^a\ki)\ 2 f(e ki )[l-f(e kf )} 
x S(ekf - eki - Hlo) (10) 



6 



where the indices i(f) count bands of initial (final) transi- 
tion states, the factor 2 takes the different spin orientations 
into account. The observed interband peaks thus allow for 
an extra check of the band structure calculations. The scat- 
tering processes smear the interband spectra, which is of- 
ten roughly accounted for by a spectral convolution with a 
Gaussian. 

Energy (eV) 



1 2 3 4 5 6 




Wavenumber (10 3 cm ') 

Fig. 7. Left panel: the anisotropic interband optical conductivity of 
MgB2 (adapted from Ref. I19p obtained by ab initio calculations. 
Right panel: the band structure close to the Fermi level (adapted 
from Ref. [3]). The origin of the most prominent interband peaks A, 
B and C is show schematically by arrows. 

Several groups [8|19l42|44j calculated a IB (u) for MgB 2 
all in reasonable agreement with each other. Fig [7] shows 
the calculated in-plane and c-axis interband conductivities 
from Ref .[19 . A remarkable anisotropy is present: there 
is no sizeable optical intensity along the c axis below 4 
eV, whereas there are two peaks at 0.3-0.5 eV (marked 
as A) and 2.4 eV (marked as B) for polarization parallel 
to the boron planes. Peak A corresponds to the transition 
between two close a bands. Peak B is due to a transition 
from the a band to the it band close to the M point, where 
a van Hove singularity strongly enhances the density of 
states. The quasi-2D character of the a bands explains why 
the transitions A and B are not optically observed along 
the c axis. The c-axis conductivity features a remarkably 
intense interband peak at about 5 eV, which comes from a 
transition between almost parallel bands. 

The a — > 7r transition (peak B) is observed in Refs.[23] 
and 25j. In the cllipsometric study by Guritanu et al. [23], 
this peak was found to be centered at 2.6 eV (FigfTJ), which 
is slightly higher than the theoretical value 2.4 eV [TI5] . The 
mismatch in the peak position suggests that the separation 
between c-band and 7r-bands is bigger than predicted by 
the theory by about 200 meV. Interestingly a shift of the 
same amount brings the results of the dHvA experiments 
[TTj in MgB2 close to theoretical predictions [IS] . Kakeshita 
et al. |25j reported the interband transition even at a higher 
energy of 2.8-2.9 eV (Fig]3]) which could be either due to the 
sample difference, or due to the fact that the peak was close 
to the experimental frequency cutoff of that study. No shift 
of the energy of this transition with carbon substitution 
was observed. 

As a result of its large optical strength and the proxim- 
ity to the screened plasma frequency, the a — > it transition 




Fig. 8. Two photographs of the same crystal of MgB2 (ac-plane) 
made in polarized light: E \\ ab (left) and E \\ c (right). Adapted 
from Guritanu et al. 1231 . 

broadens the in-plane reflectivity plasma edge and shifts it 
down by about 0.5 eV with respect to the one along the 
c-axis (as was noticed in section [3] the bare plasma fre- 
quency appears to be almost isotropic!). A beautiful man- 
ifestation of the plasma edge anisotropy is the multicolor 
appearance of magnesium diboridc [23 27]. The ac-orientcd 
single crystals looks golden in light polarized along the c- 
axis and blueish-silver for the in-plane polarization (FigJS]) . 
In the raw reflectivity data, this transition causes a second 
pseudo plasma edge (FigfTJ), that was already seen in the 
first single crystal spectra by Perucchi et al. [46] . 

Less clear is the situation with the low-lying a — > a tran- 
sition. According to the calculations, it should be seen as a 
minute dip of about 1 % on the in-plane reflectivity [53] . Di 
Castro et al. reported such a feature at about 0.5 eV and 
found that it shifts to lower energies with carbon and alu- 
minum doping. However, this effect is not obvious in other 
optical studies, for example in Rcfs. 23 25 18J . More exper- 
iments are needed to clarify this issue. 

Unfortunately, the experimental range of the c-axis opti- 
cal studies [23125] did not allow seeing the 5 eV peak (peak 
C in Fig[7]). An ultraviolet c-axis experiment is highly de- 
sirable to compare its position with theoretical predictions. 

6. Far-infrared spectra: one gap or two? 

Far-infrared transmission data of thin lead films [47] 
showed the existence of the superconducting gap before the 
advent of the BCS theory. Not surprisingly, many groups 
examined the unusual superconductivity in MgB2 with far- 
infrared spectroscopy. Although infrared spectra do not tell 
as directly as tunneling spectra about the electronic density 
of states, the superconductivity induced correlations have a 
dramatic effect on radiation absorption on the energy scale 
of the gap. In an isotropic s-wave BCS superconductor, the 
real part of the optical conductivity o\ (w) below 2A van- 
ishes at low temperatures, because photons with a lower 
energy cannot break Cooper pairs. The inductive compo- 
nent (Ji{iS) shows ~ 1/u) behavior due to the formation of 
the condensate. 

An early far-infrared grazing-incidence reflectivity 
measurement by Gorshunov et al. [48] on a porous 
polycrystalline sample of MgB2 (T c =39 K) revealed a 
superconductivity-induced rise of reflectivity with a max- 
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Table 1 

A summary of the superconducting gap values of MgB2 extracted 
from different far-infrared experiments. Values marked by * are de- 
termined from the single-gap fit, by ** - from the maximum of 
Rs/Rn, by *** - from the onset of absorption. 



sample 


rcfl/tran 


T c (K) 2A (meV) 


2A 

k R T a 


Ref. 


thin film 


tran 


33 


5.2* 


1.8 


[50] 


thin film 


tran 


30.5 




1.9 


[52] 


thin film 


refl 


32 




1.1 


[49] 


thin film 


refl 


30 


5.2* 


2.0 


M 


thin film 


refl 


35 




1.(5 


[53] 


polycrystal 


refl 


39 


3-4** 


0.9-1.2 


[48] 


polycrystal 


refl 


39 


4.8* 


1.4 


[55] 


thick film 


refl 


39.6 




1.5 


m 


mosaic 


refl 


38 


y ^fi ^ 


1.2 




single crystal 


refl 


38.5 


■i* 


1.2 





imum effect at about 3-4 meV, which was tentativeiy 
attributed to the lowest value of 2A. Further experiments 
on thin films [49|50|51|52|53"l54] . thick films [32], single 
crystals Hi 2Tj and dense polycrystals [19155] provided 
much more detailed and quantitative spectral informa- 
tion that in general supported such assignment (Tabled]). 
Note that this gap value gives an anomalously small ratio 
2A/fcsT c ~ 1 — 2 as compared to the weak-coupling BCS 
result of 3.52. However, it agrees well with the smallest 
gap observed in the tt band by tunneling spectroscopy 13J 
and ARPES [16]. 

While a gap-like feature at an energy scale of the order 
of the smaller gap A n was universally seen, no clear evi- 
dence of the larger gap (~ 7 meV according to tunneling 
[15] , ARPES [IS] and Raman [T7] data) was found in the 
far-infared experiments. This is perhaps the most puzzling 
and intriguing issue raised by the infrared spectroscopy of 
magnesium diboride. 

To simplify the comparison of different measurements, 
we begin with a basic simulation of the expected effect of 
the two-gap structure on the in-plane infrared spectra of 
MgB2 , similar to the one made by Lee et al. [51] and Lobo 
et al. [54] . We assume that a and 7r charge carriers respond 
independently to the external radiation; the condition that 
justifies using Eq.Q. Each component is characterized by 
its own in-plane plasma frequency w PiQ , scattering rate 7 
(which is mostly due to impurities at T < 40 K) and super- 
conducting gap A. We use the Mattis-Bardeen formalism 
[56] for an s-wave isotropic BCS superconductor extended 
to the case of finite impurity scattering by Zimmermann et 
al. [57] , The standard Fresnel equations are employed. 

Figure [9] shows the superconducting to normal-state ra- 
tios Js/Jn = /(0)//(T c ) for the real part of the optical 
conductivity (top), the reflectance of a bulk sample (mid- 
dle, solid line), the one of a thin film with the thickness 
of 30 nm (middle, dashed line) and for the transmission of 
the same thin film (bottom). In the first simulation (left 
panels), it is assumed that the gap values A a — 7.1 meV 



Energy (meV) 
5 10 15 20 5 10 15 20 




Wavenumber (cm 1 ) 



Fig. 9. A simulation of the far-infrared properties of a two-gap BCS 
superconductor using the Mattis-Bardeen formalism extended to ar- 
bitrary purity [57] . The two cases are considered (left vs. right pan- 
els). The values of 7 and A for a and 7r bands are given on the 
top panels. In both cases, u) PttTa =i.l eV, LO v ^ a =G.9 eV. The vertical 
lines denote 2A„- and 2 A CT . 

and A^ = 2.2 meV are the same as in the tunneling spec- 
tra [13]. The scattering rates 7^ = 12.5 meV and 7^ = 85 
meV are taken from the optical study [23] . The plasma fre- 
quencies are borrowed from the LDA calculation [5] which 
are close to the recent optical results [23125) as discussed 
in Section[5J We ignored the effective mass renormalization 
due to electron-phonon interaction which does not qual- 
itatively change the results of this simulation. The con- 
ductivity ai(u>) shows two independent thresholds at 2A n 
and 2 Act. The reflectivity ratio is slightly less than 1 for 
lo > 2Ao- showing a minimum somewhat above 2A CT , but at 
lower energies it suddenly grows and saturates (for a thin 
film) or shows a broad maximum (for a bulk sample) at 
about 2A7r. The film reflectance structure becomes quite 
sharp if the thickness is less than the penetration depth. 
This generally makes thin films experiments more sensitive 
to the gap value than the ones made on bulk samples. The 
film transmittance shows a pronounced peak centered near 
2A CT , while the n gap only slightly modifies the curve shape 
without producing any sharp spectral features. 

Jung et al. [50] measured the far-infrared transmittance 
of a thin (~ 50 nm) MgB 2 film (T c = 33 K) with the c axis 
perpendicular to the surface fFig llOp . The transmittance 
in the superconducting state showed a peak at around 6.5 
meV at 5 K. The best fit assuming a single gap value was 
achieved with 2A ps 5.2 meV. The temperature dependence 
of the gap was close to the prediction of the BCS theory. 
However, the peak position is about two times (!) lower than 
the prediction of the two-gap model based on the tunneling 



8 




Fig. 10. Far-infrared transmittance of a thin c-axis oriented MgB2 
film (T c as 33 K). Solid line is fit by a one-gap model. The inset shows 
the temperature dependence of the peak position and 2A. Adapted 
from Jung et al. 1501 . 

result (FigEl). To be more accurate, one needs to correct 
for a reduced T c by taking a smaller value of A CT . From 
the solution of full Eliashberg equations, Dolgov et al. [58] 
found that A CT should be about 5.5 meV for T c =33 K, which 
is close to the recent observations on neutron-irradiated 
MgB2 samples [S5] . However this correction seems to be by 
far insufficient to fix the mismatch. 




2 4 6 8 10 12 14 16 
Photon Energy (meV) 



Fig. 11. Real part of conductivity ci(aj) of thin MgB2 film (T c = 30.5 
K) normalized to its normal state value at 40 K. From top to bottom: 
50, 30, 27, 24, 17.5 and 6 K. Solid curves obtained by Mattis-Bardeen 
calculations (27, 24, 17.5 and 6 K). Inset: real (circles) and imaginary 
(squares) part of normal state conductivity at 40 K, along with a 
Drude calculation (lines). Adapted from Kaindl et al. |52] , 

Kaindl et al. [52] reported &\(uj) and <j 2 (uj) of two MgB 2 
films (T c = 30 and 34 K) directly obtained by terahertz time 
domain spectroscopy in transmission geometry. The trans- 
mission peak was at about 7 meV for both films, which is 
close to the observation in Ref. 50J. A strong depletion of 
absorption was observed below ~ 5 meV at 6 K, although 



the data do not exclude that the full conductivity suppres- 
sion is at somewhat lower energy (Fig fjTj) . 



R 







MgB 2 


s 




T C = 32K 

J0 




>'\ 

\ \ 




1 1 






"* 'i 

A J 1 


\# 

• • 5K 

r >\ • 20 K 

Vtf • 25K 

* 40K 



Fig. 12. Far-infrared (lines) and submillimcter (dots) reflectance of 
a thin MgB2 film. Adapted from Pimenov et al. |49| . 

Pimenov et al. [3H] obtained the reflectance of a thin 
MgB 2 film with a similar T c (32 K) (Figfl2]). The spectral 
range was extended down to 0.5 meV by using a quasiopti- 
cal backward wave oscillator technique which gives directly 
rji and <72- The onset of absorption was observed at about 
3 meV, which was ascribed to the value of 2 A. Although 
this number is smaller than 4.2 meV found in Ref. [50] . one 
should realize that the way of determination is different. In 
fact, a pronounced reflectivity dip is observed in Ref.[49j 
at ~ 7 meV, which is close to the transmittance peak from 
Refs. [50] and [S2]. A similar result was obtained by Lobo et 
a/. [54] on a c-axis oriented film with T c « 30 K. This peak- 
dip match suggests the mutual consistency of different thin 
film studies. 

While the straightforward application of the model as- 
suming a large gap separation clearly fails to describe the 
existing infrared spectra, one may wonder whether these 
data are still compatible with the presence of two, less dis- 
tinct, effective gaps. Lobo et al. [54] found that a two-gap 
model with A^ = 1.8 meV and A CT = 3.5 meV improves the 
fit of the reflectivity curve compared to the best fit by the 
one- gap model with 2 A = 2.6 meV, although the improve- 
ment was not sufficient to conclude that two gaps must ex- 
ist. In Fig[9](right panel) a second simulation with the same 
formulas is shown where the parameters are changed in two 
ways: (i) the gaps are reduced to A^ = 1.8 meV and A CT 
= 3.5 meV as in Ref. [54] . which nevertheless leaves them 
quite distinct and (ii) the sample is assumed to be far in 
the dirty limit (7^ = 125 meV and 7^ = 250 meV), which 
is perhaps a better description for thin films as compared 
to single crystals. This model qualitatively matches both 
reflection [49)54] and transmission |50|52] spectra. 

Single crystals of magnesium diboride have an advantage 
of higher T c and smaller residual resistivity as compared to 
thin films. The first single-crystal optical measurement of 
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MgB2 was made by Perucchi et al. [H] on a mosaic of ab- 
oriented samples (T c =38 K). The spectra were taken down 
to 3 meV while changing temperature down to 1.6 K and 
applying magnetic fields up to 7 T, which allowed a fully op- 
tical determination of the temperature dependence of the 
c-axis critical field H C 2 (~ 5 T at lowest 1.6 K). The authors 
[46 observed the absorption threshold at about ~ 4 meV, 
which is close to the observation on thin films 52 49 . Di 
Castro et al. [24', measured the ratio Rs/Rn on separate 
ab-oriented single crystals of MgB2 (T c =38.5 K) using the 
synchrotron radiation between 3 and 12 meV to compen- 
sate the loss of signal due to the small sample size. They 
observed a relative reflectivity increase below T c , which is 
much smaller in size as compared to the one from Ref . |46j , 
but much more close to the expectation according to the 
Mattis-Bardeen-Zimmermann calculation (FigEJ) . The rise 
of reflectivity was observed only below 5 meV at zero dop- 
ing and 8 meV for a carbon substituted sample, which is 
much lower than the prediction of the model that assumes 
a large value of A a (FigUl left panels). The reason for a 
significantly different size of the superconductivity induced 
increase of reflectance in Refs.[46] and [24] is not fully clear, 
although it may be related to mosaic effects in the first mea- 
surement or diffraction effects due to a small sample size 
in the second one. Further experiments on single crystals 
may clarify this issue. 

Let us now speculate on the reasons for the failure of the 
straightforward two gap description of the infrared spectra 
in MgB2. One possibility is that the assumption about two 
independent responses of a and 7r bands has to be aban- 
doned. Another weakness of the modeling of Fig|9]is the as- 
sumption about the local electromagnetic response and the 
validity of the Fresnel equations. This would be well justi- 
fied in the dirty limit, which is likely to be applicable to the 
7r bands, but not necessarily to the a bands. The last possi- 
bility is that there are two gaps in the infrared spectra, but 
they are effectively reduced for some reason with respect 
to the ones observed by other probes. Certainly more the- 
oretical and experimental research is needed to solve this 
paradox. 

To get the full picture of the electromagnetic response of 
MgB2 in the superconducting state one would have to com- 
plement the discussed spectroscopic data with the temper- 
ature dependence of o\ and 02 (usually expressed via the 
penetration depth A) as measured by microwave 60 61 62J 
and submillimeter [63 techniques. Interestingly, the struc- 
ture of the so-called coherence peak is highly anomalous 
for a single gap superconductor, revoking the explanation 
in terms of two gaps (some of these measurements are re- 
viewed in Ref. [53]). Optical studies of photoexcited quasi- 
particle dynamics [65154] provide another useful piece of 
information. The discussion of these latter experiments is 
unfortunately beyond the current review. 



7. Conclusions 

Many optical properties of magnesium diboride now ap- 
pear to be consistent with the ah initio electron and phonon 
band calculations. The long standing problem of an appar- 
ently small plasma frequency is solved due to the recent 
measurements on single crystals: it is close to the theoreti- 
cal prediction of 7 eV and is almost isotropic (which should 
be regarded, in fact, as a coincidence, since the plasma fre- 
quency of the a bands is very small along the c axis). The c- 
axis plasma frequency strongly decreases with carbon dop- 
ing, in contrast to the one parallel to the boron planes. The 
observed energy of the a — > tt interband transition (2.6 
eV) suggests that the LDA calculations underestimate the 
separation between a and ir bands by about 200 meV, in 
agreement with the de Haas- van Alphen experiments. Op- 
tical data in a combination with DC resistivity curves show 
that the a bands are characterized by a stronger electron- 
phonon coupling but a smaller impurity scattering as com- 
pared to the tt bands. 

Far-infrared spectra clearly show a gap-like onset of ab- 
sorption indicating 2A = 3-5 meV and an anomalously 
small ratio 2A/ksT c =l-2. This value likely refers to the 
gap in the ir bands. A serious problem which remains is 
the absence of a clear far-infrared signature of the large 
gap seen in other experiments in MgB2 . Before this issue is 
resolved, our understanding of two-gap superconductivity 
can hardly be considered as a complete. 
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